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Abstract. At large collision energy y/s and relatively low momentum transfer Q, one 
expects a new regime of Quantum Chromo-Dynamics (QCD) known as "saturation" . 
\Q ' This kinematical range is characterized by a very large occupation number for gluons 

inside hadrons and nuclei; this is the region where higher twist contributions are as large 
as the leading twist contributions incorporated in collinear factorization. In this talk, 
' I discuss the onset of and dynamics in the saturation regime, some of its experimental 

I/"*) . signatures, and its implications for the early stages of Heavy Ion Collisions. 
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1. Hadrons and nuclei at high energy 
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A nucleoli at low energy can be seen as made of three valence quarks, constantly 
interacting via gluon exchanges and virtual fluctuations at all space-time scales smaller 
than the size of the nucleon itself. In an interaction process with a probe, only those 
fluctuations which are longer lived/larger than the resolution of the probe are actually 
relevant. In addition, interaction processes at low energy are made very complicated 
by the fact that the constituents of the nucleon can interact during the time seen by 
the probe. When one boosts the nucleon to a higher energy, all its internal time scales 
are dilated, which simplifies the interaction process with the probe@: the interactions 
among the constituents of the nucleon now occur over much larger time-scales, and 
therefore the probe sees only a collection of free constituents. Moreover, the life-time 
of the quantum fluctuations is also time dilated, and thus the number of gluons taking 
part to the interaction process increases with the collision energy. Simultaneously, the 
fluctuations that were already important at the lower energy are now evolving so slowly 
that they can be considered static over the time-scale seen by the external probe. 

However, this growth with energy of the number of gluons in the wave-function 
of a nucleon (or nucleus) cannot continue indefinitely. Indeed, it would imply that 
nucleon-nucleon cross-sections grow faster than what is allowed by Froissart's unitarity 
bound. In fact, an important aspect of the physics is missing in the above picture: 
gluons can recombine when their occupation number is large, a process known as gluon 
saturation pQ. To quantify when this new phenomenon occurs, one must compare the 

* We assume that we are in a frame in which the probe has not changed 
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Figure 1. Saturation domain in the (Q 2 , x) plane. We have also represented the larger 
domain where geometrical scaling occurs. 



number of gluons per unit of transverse area^l, p ~ xG(x } Q 2 )/ttR 2 , and the cross- 
section for recombination, a ~ a s /Q 2 . Saturation occurs when 1 < per, or equivalently 
Q 2 < Q 2 s { x )i where Q s (x) = a s xG(x,Q 2 )/7iR 2 is known as the saturation momentum. 
The equation Q 2 = Q 2 (x) delineates the border of the saturation domain, indicated in 
figure [H Phenomenologically, Q 2 varies UbflA'/^ with the nucleus atomic number 
A and the momentum fraction x. 

2. Color Glass Condensate 

The Color Glass Condensate (CGC) is a description of the nature of the saturated 
nuclear matter, in which a saturation scale emerges naturally. In the CGC description, 
the fast partons (large x) are described as static color sources - represented by a density 
p{x±) in the transverse plane - that act as external sources. Conversely, the slow 
partons (low x), dominated by gluons, are described by the usual dynamical gauge fields 
N 1 . Such a separation of the degrees of freedom was first proposed in the McLerran- 
Venugopalan model [2]. In order to make this description complete, one needs to specify 
the distribution W y [p] for the hard sources in a projectile that has evolved up to 
the rapidity Y = ln(l/x). This distribution is in principle non-perturbative, but its 
evolution with rapidity Y is controlled by an evolution equation that can be derived in 
perturbation theory [3], 



known as the JIMWLK equation. This functional evolution equation can be seen as an 
extension of the BFKL equation [5J that incorporates the effect of gluon recombination. 
It has a very useful mean field approximation: the Balitsky-Kovchegov equation [I]. 
These evolution equations are presently known at leading logarithmic accuracy, but 

§ This is the relevant density at high energy. Thanks to Lorentz contraction in the longitudinal 
direction, soft gluons belonging to different nuclcons have overlapping wave- functions and act coherently. 
II The A dependence follows from the fact that all the gluons at a given impact parameter act coherently, 
and the x dependence can be inferred from the gluon distribution measured at HERA 



dW Y [p] 
dY 



n[ P ]w Y [p] 



(i) 



QCD at small x and nucleus-nucleus collisions 



3 



several recent works have extended them in order to include the Next-to-Leading Order 
corrections that come from the running of the strong coupling constant [6]. Another 
direction actively pursued is aimed at including "pomeron loops" in this description (see 
[7] for recent advances). Pomeron loops can be seen as fluctuations that are important 
in the region when the gluon occupation number is low, and may affect crucially the 
rapidity evolution of scattering amplitudes since it is controlled by the tail of the gluon 
distribution. 

3. Experimental evidence of saturation 




Figure 2. Left: geometrical scaling in the 7*-proton cross-section at HERA. Right: 
Limiting fragmentation in RHIC data. 

The first evidence of saturation comes from the "geometrical scaling" observed 
in Deep Inelastic Scattering at HERA [8]. There, the measured j*p cross-section at 
x < 0.01 - which in principle depends on both variables x and Q 2 - was found to 
depend only on the combination r = ln(Q 2 /Q 2 (x)) with Q 2 S oc x~ ' 3 (see the left panel 
of figure [2]). Such a scaling is a direct consequence of the behavior of solutions of the 
evolution equation for the distribution of sources in the target proton, and the value 
of the exponent 0.3 can be accounted for from a careful analysis of these solutions [9]. 
Going beyond this simple scaling, saturation physics leads to very good fits to a broad 
set of data from HERA and RHIC gU\. 

Another experimental result which is suggestive of saturation is the so-called 
"limiting fragmentation" (right panel of figure [2]). By shifting the rapidity axis by 
the beam rapidity, one can see that the rapidity distribution for produced particles at 
collisions of various energies tend to some universal curve in the fragmentation region 
[TT] . In the CGC framework, this property follows naturally from the unitarization of 
scattering amplitudes in the dense target, and the approximate Bjorken scaling in the 
fragmenting nucleus. The limiting curve then appears to be a reflection of the parton 
distribution at large x in the fragmenting nucleus [T2] . 

The observation of the suppression of hadron spectra at rather large p± and forward 
rapidity in deuteron-gold collisions at RHIC [13] (see figure |3] for BRAHMS results) 
provides further information on the onset of saturation. The fact that the nuclear 
modification factor R^a is larger than unity at mid rapidity is believed to be an effect 
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Figure 3. Forward suppression observed by BRAHMS. 



of multiple scatterings (Cronin effect), and the suppression of this ratio at forward 
rapidities is a consequence of the shadowing that builds up via the evolution in Y [14J. 

4. Initial particle production in nucleus-nucleus collisions 

In order to use the CGC framework to describe the collision of two nucleons or nuclei, one 
needs two color sources px^, representing respectively the two projectiles, that couple 
to the gauge fields (see [15] for a detailed discussion on this) via the following current, 

J» = 8»+5(x-) Pl (x ± ) + 5^5(x+)p 2 (x ± ) . (2) 

As illustrated in figure HI the description of Heavy Ion Collisions is usually split 
into several stages, each of them being described by different theoretical tools. It 
is expected that perturbation theory is appropriate for the very early stages, during 
which the constituents of the incoming nuclei are released and particles are initially 
produced. However, in a typical central heavy ion collision at high energy, 99% of the 




Figure 4. Left: the successive stages of a heavy ion collision. Right: typical diagram 
contributing to gluon production in the collision of two nuclei at high energy. 

produced particles have a transverse momentum below 2 GeV - i.e. below the saturation 
momentum. In this region, collinear factorization breaks down, and it is believed that 
the CGC framework is more appropriate for describing the initial particle production. 

In the saturation regime, the sources p\p are typically of order of 1/g, which alters 
the power counting for gluon production: adding an extra source to a given diagram 
does not change its order, implying that the color sources px,2 should be included to 
all orders. A typical diagram contributing to gluon production from these sources is 
illustrated in the right panel of figure HI In general, there are several disconnected 
subdiagrams, including some that do not produce any gluon (vacuum diagrams) [PT] . 
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Figure 5. Gluon multiplicity at leading order. Left: tree retarded diagrams 
contributing to the gluon spectrum. The color sources are represented by red and 
green dots. Right: reformulation as an initial value problem on the light-cone. 

At leading order [16], one must sum all the tree diagrams, with an arbitrary number 
of color sources - in general a very complicated task. However, for the inclusive gluon 
spectrum, this can be done from the classical solution of Yang-Mills equations with 
retarded boundary condition^. The retarded nature of the boundary conditions leads 
to an important simplification: one can reformulate the problem as an initial value 
problem where one specifies the fields and their first time derivatives just above the 
light-cone (and these initial conditions can be obtained analytically). This classical 
solution of the equation of motion is a sum of tree diagrams with retarded propagators, 
as illustrated in figure El An important feature of the classical solution of Yang-Mills 
equations is its invariance under boosts in the longitudinal solution, which implies that 
the field has only zero modes with respect to the rapidity variable rj. 

4-2. Gluon production at Next to Leading Order 

In order to make the foundations of this description more robust, it is important to 
study loop correctionsEl One reason for doing so is to check the factorization of the 
leading logarithms of 1/x in the evolution of the distribution of sources Wfpi^], which 
is necessary for the internal consistency of this framework. Another motivation comes 
from the recent observation that the LO boost invariant solution is unstable against 
rapidity dependent perturbations [19] - and loop corrections are a natural origin for 
such perturbations. 

Although the LO calculation of the gluon spectrum is non-perturbative since it 
requires to solve the Yang-Mills equation to all orders in the sources pi^, it turns out 
that one can express the 1-loop contributions as a perturbation of the initial value 
problem encountered at LO. There are two contributions, illustrated in figure [6j The 
NLO correction can be obtained by the action of some analytically calculable operator 

p Less inclusive observables can be expressed at LO in terms of classical solutions with more complicated 
boundary conditions [17] . 

+ Quark production - also a 1-loop contribution - has been evaluated in [18) . 
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on the LO contribution, seen as a functional of the fields and canonical momenta on the 
light-cone [20] . 



SN A 



5A(x)T x + 



1 



X(x,y)T x T y N LO [A in [ Pl ,p 2 ]] . (3) 



ccglight— cone £c,j/Glight— cone 

In this formula, 5A(x) and Y,(x,y) are respectively the 1-point and 2-point functions 
represented in green in figure [6] - they are in principle calculable analytically. T x is the 
generator of shifts of the initial condition at point x on the light-cone for the solution 



of Yang-Mills equations^, i.e. a functional derivative S/SAi n (x). Eq. ([3]) is sufficient to 





Figure 6. The two contributions to the gluon multiplicity at NLO. 

discuss the structure of singularities that arise at 1-loop : 

i. the coefficients 5A and S are infinite, due to divergences in integrations over 
rapidity [2T] . 

ii. the quantity T x A(t, y) diverges when the time r goes to infinity. This is due to 
the instability of the boost invariant classical solution found in [TP] . 

The first type of divergences is related to the evolution with rapidity of the 
distributions W^pi^]- For the CGC framework to be consistent, one should be able to 
absorb these divergences in the rapidity dependence of W^fpi^] - a statement very similar 
in spirit to the factorization of collinear divergences in conventional perturbative QCD. 
Another way to state the same problem is to set fictitious boundaries in rapidity, Yq and 
Yq, between the observable (here the number operator) and the color sources. Now, all 
integrations over rapidity in the observable are over the finite range [Yq, Yq] , and thus lead 
to a finite result. Thus, 8 A and T, become finite, but Yq and Y ' dependent. Naturally, the 
distributions of color sources, PFfpi^], must now be evolved from the beam rapidities 
to Y and Y ' respectively. The factorization alluded to before is equivalent to the 
requirement that the Yq and Y ' dependence cancel between the 1-loop correction to the 
observable and the rapidity evolution of the distributions of sources. For this cancellation 
to happen, one will need to establish the following relation for these divergent terms 
EG] : 



6N K 



divergent 
coefficients 



(Y -Y)rt M + (Y - Ytftf [p 2 ] N LO [A 



(4) 



* Eq. ([3]) is only a sketch of the actual formula. Indeed, N LO depends on both the initial fields A- u 
and momenta IIj n on the light-cone, and there should therefore be shift operators for both of them. 
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where TC[p] is the same Hamiltonian as in the JIMWLK evolution equation. 

The second type of problem is due to an instability of the boost invariant solution 
of the classical Yang-Mills equations |19j . After resummation of these divergences to all 
orders, the operator acting on N LO in Eq. (J3J) will be replaced by some functional Z[T X ]. 
Formally, one can perform a "Fourier transform" of this functional by introducing some 
auxiliary field a(x)on the light-cone§, and then rewrite the contribution of the resummed 
unstable modes as§ 

J [Da{x)\ Z[a] exp {z J a{x)T x } N LO [A in ] 

= j [Da(x)\ Z[a] N LO [A m + a] . (5) 

In other words, resumming the unstable modes can be done by solving the classical 
Yang-Mills equations with a fluctuation a(x) superimposed to the boost invariant 
initial condition A- in) with a distribution Z[a\. This is the formal justification for the 
fluctuations added to the initial conditions in [Tjj]. By a completely different approach 
|22j, it appears that the functional Z[a] is a Gaussian, suggesting that the resummation 
of the unstable modes is simply an exponentiation. 



SN NLO 



unstable 
modes 



5. Conclusions 



QCD at small x, and in particular gluon saturation, has become an important aspect of 
hadronic collisions at RHIC and LHC energies. In nucleus- nucleus collisions, the CGC 
is the appropriate framework to study the early stages of the collision, because most 
of the particles are produced from low x gluons. It seems now technically feasible to 
compute the 1-loop corrections to the gluon spectrum, and to resum their diverging 
terms. Doing so would allow one to prove for the CGC framework factorization results 
that are crucial for the overall consistency of this approach. 

Acknowledgements : Important parts of this talk are based on work done in 
collaboration with K. Fukushima, T. Lappi, L. McLerran and R. Venugopalan. 
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